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SUMMARY 


An analytical and experimental investigation was made to 
determine the possibility of predicting the static and rotary 
directional stability derivatives of a flying boat of the con- 
ventional short-afterbody type. Starting with the angle of trim 
obtained from a tank test at a specified load and speed, the ex- 
ten” °* the wetted areas of forebody and afterbody is computed 
on the basis of available empirical relations available for 
planing processes. It is shown that for a range of pre-hump 
speed the lateral forces due to the yaw can be obtained by con- 
sideration of the inequality of loading on two sides of the bot- 
tom due to the change of the effective angles of attack resulting 
from vectorial addition of the angles of deadrise, yaw, and trim. 
It is also shown that this action results in identical position 
of centers of pressure for the vertical and lateral forces. Qn 
this basis, the static and rotary stability derivatives for yaw- 
ing motion are computed. 

The experimental verification of the method in regards to 
static derivatives was made for three operating speeds at Cy 
of 2.63, 2.6, and 2.34. A reasonable agreement of the test and 
analysis data is shown for the two higher speeds, but the agree- 
ment at the lowest speed is uncertain, probably because of the 
lack of necessary data on the configuration of the wake of the 
forebody. Such data are presently available only for Cy Gn 


3.0 and above, 





The work was done at the Graduate School and Experimental 
Towing Tank of Stevens Institute of Technology in partial ful- 


fillment of the requirements for the degree of Master of Science. 





INTRODUCTION 


Locke and Pierson (References 1 and 2) among others point 
to the ever-increasing instances of flying-boat directional in- 
stability in the past years with an increase in operating loads, 
pe oeei ally in the pre-hump region. Many investigations on spe- 
cific types of hulls have been carried out in order to determine 
what kind of modifications to the hull will give the best direc- 
tional stability. Extended afterbody and planing-tail hulls 
(Reference 3) have received quite a bit of attention. Besides 
considering the type of hull which will best give the desired di- 
rectional stability, much of the present work has been directed 
toward auxiliary devices which will aid directional control in 
the pre-hump region. Such devices include water-rudders (eee 
ence |), reversible-pitch propellers and hydro-flaps (Reference 
5). Hydro-flaps appear as an excellent means of obtaining good 
directional control. From current trends it is entirely possible 
that the problem of directional control in the pre-hump region 
will be neatly by-passed in favor of auxiliary control devices 
and little attention will be paid to the basic hull unless the 
hull should exhibit some extremely radical directional tendencies. 

It is the intent of the present work to return to the 
basic hull without any auxiliary control devices and to determine 
what the hull stability will be from data obtained in ordinary 
planing test runs. The inherent directional instability of the 


forebody is due to the excessive wetted length extending forward 





of the C.G. of the forebody (Reference 6). Therefore the after- 
body must supply the necessary stabilizing moment for the entire 
hull to be directionally stable. The contribution to stability 
of the afterbody is extremely difficult to predict because of the 
changes in the wave form in the wake brought about by changes in 
trim and speed of the forebody, The yawing moments due to the 
afterbody directly depend, of course, on this shape of the wake. 
The moments created by the afterbody can be considered to be of 
two types. The first is due to the shape of the afterbody and the 
shape of the wake while the second is of an accidental nature. 
That is, moments of the second type are due to water flowing over 
the curved sides of the afterbody or over various portions of the 
tail cone. Attention is concentrated only on the first type 
since the second type is not amenable to any simple analysis, 
being dealt with in practice by the use of stringers and other 
type of flow-corrective devices. 

There is quite a bit of available information on porpois- 
ing and longitudinal stability, but very little information con- 
cerning directional stability, especially from an analytical 
viewpoint. This is due to the complex nature of the flow exist- 
ing around the afterbody as already suggested. In the early 
stages of the present investigation, it was felt that the range 
of investigation should be limited to the pre-hump region where 
most of the instability of the seaplanes is encountered. In the 
region between hump and getaway, the hull has either sufficient 
inherent directional stability or aerodynamic and thrust control 


sufficient to maintain a straight course. The range of investi- 





gation was also limited to the horizontal plane. That is, heel 
angle, rolling, and pitching velocities were neglected. 

The general problem of the pre-hump region involves oth 
low-speed taxiing as well as pre-hump speeds involved during 
take-off. This report as well as being limited to the pre-hump 
region also excludes low-speed taxiing. This choice has been 
made since during pre-hump speeds used for take-off, the sides 
are essentially dry and most of the weight is supported by dy- 
namic loads on the forebody and afterbody and therefore use can 
be made of available planing theory. In the case of low-speed 
taxiing, it is not possible te use the simpler planing theory, 
since the curved portion of the forebody and sides are wet. Any 
analysis for low-speed taxiing is complicated and involved. 

Using a composite chart of directional stability for a 
conventional short-afterbody hull of well-known design (Reference 
1), three speed coefficients (2.34, 2.46, and 2.63) were chosen 
for which the directional stability derivatives were to be pre- 
dicted. The static directional stability derivatives in nondi- 
mensicnal form were available for only one of the speeds chosen 
from Reference 7. Moment due to yaw and transverse force due to 
zero angle of yaw (i.e., on-course motion) were reported in this 
case for the speed coefficient of 2.62. The fundamental problem 
is to find a mathematical analysis which predicts the existence 
of certain forces on the hull, and tc show that by the existence 


of these forces the static stability already known must result. 





For the speed coefficient of 2.62, the analysis will pro- 
ceed on the following basis. The centers of pressure of both the 
forebody and the afterbody will be determined. The assumpticn 
is then made that for small deviations from a straight course 
the centers of pressure and the wetted lengths of both the fore- 
body and the afterbody do not change but remain fixed. Since 
the problem has already been confined to speed ranges where the 
sides are not wetted, and therefore since only the bottom areas 
are wetted, it may be concluded that the lateral force produced 
due to an inequality of pressures on the two sides of the hull 
acts at the same point as the vertical force produced by planing 
action. That sh the centers of pressure for lateral loads and 
for vertical loads are coincident. For the case of zero yaw 
angle, the pressures acting on either side of a vee-bottom hull 
have a vertical and horizontal component of which the horizontal 
forces, being equal and opposite, cancel each other. 

When the hull is given some arbitrary small yaw angle, a 
change in the effective angle of trim resuits which is positive 
for one side of the hull and negative for the other. Therefore 
the lifting force will be increased on one side of the hull bot- 
tom and decreased on the other. The horizontal components of 
these changing lift forces are opposite but no longer equal, and 
a resultant horizontal force will be produced. This type of ana- 
lysis is also applicable to the afterbody but with the additional 


consideration of the change of the water surface in the wake. 





AS the hull is placed at some angle of yaw, the afterbody is 
displaced and the wake is no longer symmetrical about it. The 
mean lateral inclination of the water surface in the wake will 
cause the additional inclination of the resultant force, there- 
by producing an additional lateral force. 

The sum of the horizontal forces on the forebody and 
afterbody will indicate the magnitude of transverse force pro- 
duced by a given yaw angle. The horizontal forces multiplied by 
their respective arms, as determined from the centers of pres- 
sure of the vertical planing forces produced, will give the 
magnitude of the yawing moment produced for the same given yaw 
angle. It will be shown that satisfactory agreement with known 
static stability derivatives will be obtained by this method. 

The distances from the centers of pressure of the fore- 
body and afterbody having been determined, certain relations may 
be formalized which make use of these distances and the static 
stability derivatives, in order to obtain the rotary derivatives. 
In this connection, use is made of practices common in the solu- 
tion of aircraft dynamic stability. An angular velocity creates 
a small change in the local yaw angle of the forebody and after- 
body which will in turn produce additional transverse forces. 
The transverse velocities produced by the angular velocity will 
be proportional to the distance the force is located from the 
center of gravity. The change in yaw angle at each center of 


pressure will then be a function of the ratio of the local trans- 





verse velocity and the incoming velocity. By following these 
principles, the expressions for the rotary derivatives using 
static stability derivatives may be derived as shown in the sec- 
tion on mathematical development. 

Once all the derivatives have been established, the com- 
plete dynamic stability can be determined by resorting to the 
solution of the equations of motion as developed in Appendix I. 
The dynamic stability will be seen to depend on static transverse 
forces, static yawing moments, rotary transverse forces, rotary 
yawing moments, and the hull inertia parameters. The deriva- 
tion in Appendix I is presented only in the interest of the com- 
plete picture of hull stability and no attempt is made for a 
solution of the roots. Suffice to say that positive roots in- 
dicate amplification of the disturbance while negative roots in- 
dicate a decay of the disturbance leading to stability (Reference 
e)). 

Static stability has been defined as the tendency of a 
disturbance to disappear with time with the center of gravity 
constrained to move in a straight line. Therefore, the tests 
conducted in this investigation produce a measure of static sta- 
bility since the C.G. is constrained. Under such conditions, 


four types of moment vs. yaw curves may result (Reference 9): 


Type of Stability Slope of ON'/OB 


positive stability negative 

neutral very small positive or zero 
negative stability positive 

"hooking" instability curve discontinuous at small 


angles 





Of the types of moment curves possible, the statically 
stable type were obtained for the three cases tested in this in- 
vestigation. However, the lowest speed coefficient used came 
quite close to the beginning of the region of "hooking" insta- 
bility. It is not considered unfortunate that three cases were 
chosen in which the hull used was known to be stable. The pri- 
mary intention of this investigation is to illustrate that the 
analytical method presented has application in the range chosen. 

It has been suggested in other reports (Reference 7) that 
static stability is not a good criterion of dynamic stability. 
In fact, several instances of static instability in model test- 
ing have been shown by calculation to be dynamically stable. 
Criticism against static stability suggests that it is an arti- 
ficial concept and in effect presupposes a lateral constraint 
which is not present and artificially removes from the free body 
one degree of freedom. In the present work, only the direction- 
al stability is considered, but the reader is reminded that the 
rotary derivatives are also necessary if the maneuvering charac- 
teristics are to be investigated. 

The tank tests in this investigation were conducted in 
Tank No. 1 of the Experimental Towing Tank, Stevens Institute of 
Technology, Hoboken, New Jersey. 

Acknowledgement is made to the personnel of the Experi- 
mental Towing Tank who helped in many ways with the work at all 
stages. Particular gratitude is due Professor B.V. Korvin- 


Kroukovsky for his very valuable help and patient advice. 





SYMBOLS 


Nondimensional Coefficients 
eee 


y! 


r' 


Re 


Load coefficient 

Speed coefficient 

Hydrodynamic trimming moment coefficient 
Hydrodynamic yawing moment coefficient 


Hydrodynamic transverse force coeffi- 
cient | 


Dimensionless angular velocity about 
Z-axis 


Resistance coefficient 


Dimensionless hydrodynamic force coeffi- 
cient at C.P. of forebody 


Dimensionless hydrodynamic force coeffi- 
cient at C.P. of afterbody 


Coefficient of lift (vee planing surface) 


Coefficient of lift (flat planing surface) 


Coefficient of friction drag 


Ratio of mean wetted length to maximum 
beam 


Reynolds number 


Forces and Moments 


N 


Y! 


yr! 


{] 


i 


If 


A/wb? 
V/ Veb- 
m/w! 


N/V 


v/v" 


rb/V 


R/Si*D* 
2 2 
F,/SV b 
a Wee 
F,/SV"b 


2 
20, /Cy 


Hydrodynamic yawing moment about Z-axis (moments positive 


as per right-hand rule) 


Hydrodynamic force along Y-axis (force positive to star- 


board) 


Hull resistance 


Hydrodynamic transverse force at C.P. of forebody 


Hydrodynamic transverse force at C.P. of afterbody 


= Loe 





Derivatives 

ar Static hydrodynamic yawing moment coeffi- By ON'/O 8B 
cient derivative p 

Y, Static hydrodynamic transverse force co- Y', = dY'/d B 
efficient derivative P 

N',, Rotary hydrodynamic yawing moment coeffi- N'. = ON /r' 
cient derivative 7 

Y! Rotary hydrodynamic transverse force co- Y' =o Y'/gr' 

i ae : : ie 
efficient derivative 

Other Symbols 

b Maximum beam, ft. 

V Velocity, ft./sec. 

B Yaw angle, deg.; angle between the velocity vector at C.G. 
and the X-axis (yaws positive about Z-axis in accordance 
with right-hand rule) 

iC Trim angle, deg.; angle between the undisturbed water level 
and the keel reference line (trim positive when keel line 
makes bow end above undisturbed water level) 

A Displacement weight, lb. 

R Distance from C.G. to center of pressure; subscript 1 for 


forebody, subscript 2 for afterbody 
Mean wetted length 


Chine wetted length 


d, Keel wetted length 
m Subscript for all data applying to model 
Constants 
g Acceleration of gravity (32.2 ft,/sec.*) 
W Specific weight of water (62.3 lb. /ft.?) 
Density of water (1.937 slugs/ft.-) 
v Kinematic viscosity of fresh water (1.08 x 107? ft.°/sec. 


at 20° centigrade) 


=< ee 





ANALYTICAL DEVELOPMENT AND MATHEMATICAL ANALYSIS 


The forces acting on a hull in two-step planing due to 
some given angle of yaw may be represented as shown in Figure 
l(a). That is, at some angle of yaw there will be a moment about 
the C.G. and a transverse force which may be imagined to act at. 
the C.G. 

The transverse forces and moments produced by varying the 
angle of yaw represent the static stability derivatives ON'/OB 
and OY'/AB. Actually, the true picture of the situation is 
presented in Figure 1(b), where the transverse force is seen to 
be the addition of two transverse forces, one acting at the fore- 
body center of pressure and the other at the afterbody center of 
pressure. The moment about the C.G. is the result of the fore- 
body and afterbody forces acting at their respective arms. Yaw 
angle, hydrodynamic yawing moment and transverse forces are all 
assumed to be positive. Using the notation shown, the two fun- 
damental relations will be: 

Y =F, +F (i) 


aD 
eae (2) 


u 2 


Solving these equations, an expression for the forces By and 


P, may be derived as follows: 


Eliminating F, in equations (1) and (2) gives: 


1 
Ae iie- Fa)e- Raul, 


N 


=e, + + ae 


ee as 





Substituting the value of F, in equation (1) gives 


2 
iL, Y-N 
a “Goat 
Seo 
Therefore, 
N+ hyy 


Oh ese 


- 13 8 
Therefore, ) 
eae (3) 
C ) + UF ; ) 
eae | 


| 
Forebody Planing Data 

Having obtained expressions for the forces BS and By ‘ 
it is necessary to determine the centers of pressure of the fore- 
body and afterbody. Leaving the equations for FL and Fy, ' 
previous data are analyzed for hull conditions which led to known 
values of the static derivatives. Table I gives the necessary 
model dimensions while Condition 1 in Table IV gives the basic 
knownquantities for the hull to be analyzed. 

It was necessary to begin with an assumption concerning 
the amount of the total load carried by the forebody alone, since 
the load coefficient for Condition 1 was that for the entire hull. 
Once the assumption has been made, the hull analyzed, and all 
conditions found to be in agreement, the assumption is no longer 
an assumption but becomes a fact leading to a unique solution. 

In view of the relative loads in two-step planing above the hump, 
it seemed reasonable to assume that under present conditions the 


forebody carried 81% of the total load. Therefore 


CA = »61(.8) = 648. From the equation of lift coefficient due 


— 





to an angle of deadrise, 


2Cn 
2 eee “ 
Cr 2 6.87 » 1885 e (5) 


Unfortunately, the relations involving lift coefficient and the 
aspect ratio of the wetted area are such which require trial and 
error solution if recourse is not made to nomographs available 
in Reference 6. It was found that the equations gave more satis- 
factory results because of the greater accuracy. By trial and 
error, then, the value of the lift coefficient for a zero dead- 
rise planing surface is obtained from the lift coefficient of a 


vee-planing surface: 


- ~6 
c, =a - -0065(B) C, (6) 
B O O 
6 
Cegpoeee = .,13-C. ~ 
Le Ly L, 
= ,251 - .143 (.436) 
= .251 - .0623 
= ,1887 ; 
Therefore, 
Cc, = ou 
O 


To determine the aspect ratio, the following relation is used: 


aa 2 
© = ,0120(A)2 + .0095 5 SC Cs (7) 
piel Cy 


which, for the values required, reduces to 


> 


Bk 
2 





a 
2 


O1SE = .0120(r)2 + .00138(r)* 


-~ Lh - 





Once again by trial and error, the following is obtained: 


si 


can Z 
til 


i] 


. 
.0120(1.25)2 + .00138(1.25) 
20134 + .00215 


=“EOl>>> 


fl 


Therefore, the aspect ratio becomes 
e r= eS 
Savitsky (Reference 10) has presented the following formu- 
lation for determining the position of the center of pressure of 


a vee-planing surface; 


cS = 715 = 2 eee (8) 
P C 2 
“Vy 
3.06 ee ea 
Ve 
where ; 
P 
Coyeeits 
P De 


P being the distance on the-forebody from the step. 


For the forebody being considered here, results indicate 








wnat 
C Bs he al s > aa 
: 3.06 2.00) 5 ¢ 22 
(1.25) 
elo = 15.05 + 202 
= ss oe BOS 7 
Therefore 
C, = .693 


R 


. 15 = 





By definition, 


him 


art bh (9) 
and therefore, 


- tat = 1001.25) = 19.5 £4. 


That is, the mean wetted length of the forebody is 12.5 ft. And 
from the definition of the center of pressure coefficient, 

mie 093(12.5) = 8766 Pte 
The point of center of pressure of the forebody is 8.66 ft. for- 
ward of the step. However, the dimension desired is Rs which is 
defined as the distance of the center of pressure from the C.G. 
For this investigation, the 0.G. was considered to be 3.67 ft, 
forward of the step. Therefore 


=) GOO) = sae! 


<> 
! 


ane ieclneer 


Longitudinal Wave Profile 

In order to determine the hydrodynamic force on the after- 
body, it is necessary to plot the afterbody to scale, determine 
the longitudinal wave profile, and from this determine the angie 
of trim of the afterbody and the beam. This is a critical calcu- 
lation since the load computed by reference to the longitudinal 
wave profile must balance the load supported by the forebody and 
check for zero moment about the C.G. With the aid of the after- 
body dimensions listed in Table II, Figure 9 is plotted. Korvin- 


Kroukovsky, Savitsky, and Lehman (Referencell) have compiled ex- 


= 1558 





perimental data and derived an empirical formula for plotting 
the longitudinal wave profile behind a vee-planing surface using 
the speed coefficient, angle of trim of the forebody, and aspect 
ratio as the entering parameters. Reference 1] states that the 
range of applicability of equation (10) should be restricted to 
values of the speed coefficient of not less than 3.0. It was as- 
sumed, however, that a speed coefficient of 2.62 was close enough 
to allow the use of the formulas. The following are the formulas 
eed : 

Hiaetioe+ 1 


O x 


H= -0208 “C*94 (4) (x) (10) 


H= ove] 5“ wi “ ~ ,00kk8 hme —_ 
V V 

where H is the height above the keel baseline extended in beams 

and X is the distance aft of the step in beams. Table III 

shows the solution of the formula giving the height H also in 

inches to facilitate plotting. 

It becomes necessary to introduce another correction. In 
the work at the Experimental Towing Tank, it was found that to 
make results predicted on the basis of wake survey agree with 
test results on the hull models, the wave height on the afterbody 
must be reduced to about 85%. This correction makes predicted 
wetted afterbody area agree with photographs of wetted afterbody 
areas obtained during actual test runs. After several trial and 


error computations, it was found that for the present case a wave 


2 ae 





“4 ry 
ae 5 = 
= 


height reduction to 90% gave consistent results. In the region 
ef intersection of the wave profile curve with afterbody bottom, 
the curve can be approximated by a straight line determined by 
the ordinates at the points (1) and (2) located at one and two } 
beams aft of the step. In meteine these points, the ordinates 
were reduced to 90% of the value given in Tahle IIT. The angle ; 
measured between this line and the keel determines the angle of 
trim of the afterbody. 
Afterbody Planing Data 

Figure 5 shows the dimensions necessary for determining 
the angle of sweepback of the stagnation line defining the lead- 
ing edge of the wetted area. 


This figure is based on Wagner's 


expanding plate theory, following References 6 and 10. The dis- 


tance L. has the value 
bi earl 
oe 
i ean c 


which leads to the sweepback angle defined as: | 





b/ 
tan & = aff 
all 
Therefore 
, 
+ a) ss 
caro D tans 
mn tant 


Knowing both the deadrise angle and the angle of trim, the re- 


quired angle is found to be 
tan 3 = 295 
C= 26 20" 


The angle % is 


then laid out on Figure 9 starting at the leading 





edge of the wetted length. The intersection of the constructed 
stagnation line with the chine line defines the beam of the after- 
body used in all calculations requiring the beam of the afterboay. 


The length of the keel wetted length may be measured directly as 


R= 190noen. = HONS) fh. 


and by geometry the mean wetted length equals 


Q - Ax, 8.27 ft. 
m 2 


Again the aspect ratio of the afterbody may be computed as 


— Me _ 8.27 
b 5.83 : 


The speed coefficient is changed from that of the forebody by 
the inverse ratio of the square roots of the forebody and after- 
body beam lengths, so that ; 

(Oy), = (Cydp {255 )” = 2621.52) = 543 
where the beam value constructed above equals 

bee /Ohdwee— 5,03 Tt. 
Therefore, the aspect ratio 


_ Nas Gece 


= = 9 


ae aka rm ete 
Knowing the angle of trim and aspect ratio, it is now pos- 
sible to apply equations (6) and (7) to determine the lift coef- 
ficient for the afterbody. The lift coefficient for a flat plan- 
ing surface is computed to be 
Ce eos |.0120(1.192) + .0095(.172) 
“se .0143 + 001625 


a 





_ 
= 





Cae C >? ense3) 


= 191 


while the lift coefficient fcr a vee-planing surface is shown to 


be 
Cy = ,191 - ,182(.370) 
B 
= ,191 - .068 
Saale 3 


Therefore, with the lift coefficient as computed, the load sup- 


ported by the afterbody is 


lt 


(A), = -123 Gotu EE (W7)* (5.83)° 


Gyo oO) 4b. 


Il 


Using the load coefficient, the load supported by the forebody 


is added to the afterbody load to determine how closely it agrees 


with the load coefficient of the entire hull as specified at the 


Sart . 
7 = ,048(62.3) x 104 
= 10, 300 
(A) 4 okay = 402300 + 8,950 = 49,250 
Z \ oe 
CD) pecs TO Come cceIO" = sone @ 


In order to determine the center of pressure of the after-~- 


body, equation (8) is again applied and the value of C, com-~ 





puted as 
Cp = 275 b yo 
(1.42)?/ 
- 1 é 
= 075 = 21535 + 2.2 KG » 708 ° 





In the case of the forebody, P was measured directly forward 
from the step. In the case of the afterbody, Reference 10 sug- 
gests that the distance P be measured as shown by Figure 6. 

Eee Oomeme,) = 5.06 ft. =) 70.3 ane 
Hrom Figure 9, the distance L is 127.5 in., making L/2 = 63.75 
in. The center of pressure is therefore 70.3 plus 63.75 or 
134.05 in. from the end of the afterbody. Since the length of 
the afterbody is 328.5 in., the center of pressure acts at 328.5 
minus 143.05 or 194.45 in. behind the step. The distance defined 


as Kes is then 


ee = 16.2 ft. + 3.67 ft. = 19.87 ft. 


Moment Due to Resistance and Total Pitching Moment 

The next step in analyzing the hull is to determine the 
pitching moment due to resistance. Resistance is computed for 
the model but plotted to full-scale values for convenience. Be- 
ginning with the Reynolds number which is low, the Prandtl-Karman 
relation is applied and the friction coefficient obtained. The 
friction coefficient so obtained is referred to wetted area and 
must be converted to the beam squared basis as shown. With the 
coefficient so obtained, it is then possible to determine the 
resistance in pounds. The same procedure is appiied to the 
afterbody. 


Forebody Frictional Resistance 


a ee eee 


~ 10_(12.5/22) _ ¢ 96 x 10° (11) 
1.08 x 10° 


aot ae 


_ Fey o0erie. 
(Ca), = aoe Soe sCOs31 (12) 


(C2), = Cpr sec p = 00531(2.25) (-sB5) 


00716 (13) 


Mi 


_n Py - 
R = Ch 5 Vb 1630 lb. 


Afterbody Frictional Resistance 


(Re) = 20(8:27/22) - 3.148 x 10° 
IO Bes 10 
(C_) = aot = ,00576 


) 
Fm — (3,48)*/? x 10 


_ 1 
(C.) = ,00576(1.h2) —e5 
= .00927 


R = ,00927 (2-937 | Gin aieeee) 
= 672 aioe 


The total frictional resistance is then the sum of the 
forebody and the afterbody resistance. In the process of de- 
termining the forces on the hull, it was found that the value 
of resistance was too low and therefore all values of resis- 
tance were increased by 25%. This in effect accounts for 
roughness and other unknown variations. When the resistance is 
increased by 25%, the result is found to agree very well with 
the values of resistance shown in Reference 9. Therefore the 
total frictional resistance is found to be 


1636 672 
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(Qe = 17252302) =" 2880 
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With all the forces computed for the entire hull for the 

condition being considered, it is now possible to compute the 

momenta about the C.G. of the hull. Forces considered are taken 

perpendicular to the keel reference line as shown in Figure 7. 

The inclination of the afterbody force to the vertical is con- 

sidered small and its horizontal component is therefore neglected. 


The sum of the moments yieids 


41,300(4.99) = 2880(10) + 8,950(19.87) 
206,000 = 28,800 + 177,800 
206,000 .= 206,600 


Therefore the sum of the moments about the C.G. does indicate the 
initial condition of zero moment about the C.G. to a very good 
approximation. The satisfactory agreement of both summations of 
forces and of moments is now taken as a confirmation that the 


computed 1s and ys are correct. 


Unbalance of Loads on Two Sides of the Hull 

In order to arrive at an analytical expression or proced- 
ure for the determination of the transverse forces on a hull in 
the yawed condition, it is necessary to associate a change in 
the effective angle of age on each side of the hull with a given 
angle of yaw. From Figure 8, the vee-bottom planing surface is 
considered to be equivalent to a wing with dihedral. When the 


planing surface is given a positive yaw, the force diagram shown 


for the starboard side indicates that there will be a decrease 
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in the angle of trim acting on that side, the magnitude of which 
will be 


VAa = V sin ¥Y sin 8 


Aa = Y sin 8 (for small yaw angles) . (14) 
There will be an increase in the trim angle on the port side for 
the same positive yaw angle of the same magnitude. 

Some expression must be found that will allow the increase 
in lift to be computed if the change in the angle of trim is 
known. Equations (6) and (7) may be changed slightly in view of 
the assumption that for small yaw angles the center of pressure 
and the aspect ratio remain constant. The expression for the 


original trim is then altered to include a small change in trin. 
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The original lift is defined as (Cha) 4 and the new lift 


due to an increase in the trim as (Cy ae 
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By the binomial expansion 
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Neglecting the third term as being of higher order, 
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Again the third term of the expansion is neglected as being of 


higher order. Therefore, 


(1), acumen -L) es AC =k.(k,t 7 .)°o las (1.1)(.6)-45 | 


If the value for (C, is substituted, 
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Multiplying through by k to obtain the loads gives 
eal: Sa lee 
A = A, +k — pe(neyt - (1.1) (.6)k,(k, C ) : ; 


The expression in the brackets is, to a very close approximation, 





the value of (C, > since for relatively small values of trim 
the negative term in the brackets is quite small and not serious- 


ly affected by the factor of 0.6. Therefore the complete ex- 


pression: 
A= A, +d, at 
i 
b= A fat at | (15) 


The * sign in equation (15) is put in so as to include increases 
or decreases in trim since the same result would have been ob- 
tained had the derivation been considered for a decrease in ane 
trim angle. Equation (15) states that new load on a planing 
surface is found by multiplying the original load by the ratio 
of the new trim angle to the old trim angle. 
Lateral Force Computations 

Equations (3) and (4) gave expressions for determining 
the value of the forces Fy and By acting at the forebody and 
afterbody center of pressures. The hull for Condition 1 as 
specified in Table IV has the following value for the static di- 
rectional stability derivatives: 

ee oo 


B slopes per radian ; 
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Under the assumption that the wetted length remains constant for 
small yaw angles and that the static derivativesreported are con- 
stant or linear out to 1° yaw angle, the value of the forces may 


be computed on the basis that the angle of yaw is ies 
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7 499 ft. Xt = QB = 00k b = 10 ft. 
2 = 19.87 ft. NI = NI4B = ~.003 N'b = -.03 
+ h, = 2.86 ft. V = l7 ft./sec. Aye = .0795 
Ry = .02 
t I 
Ft = =r eee cot one 
i ie 2 f 2). 86 
F, = 1.99 x 10 : pee | (u'7)*(10)* 
F, = 42h lb 
Y! = N'b 
lig 02 + .03 => 
t= = ~ = 
F, ar [ 2.01 x 10 
F, = 2.01 x 10 ae (7) *(10)* 
a lier ke 


Therefore, for the condition specified, the static deriva- 


tives will have the value shown if the forces computed act at the 
centers of pressure at the distances shown when the hull is 


placed at an angle of yaw of flee 


With the forces known that must act in order to produce 


a given static stability, a correlation must be made and shown 

to be valid if the static stability is to be predicted from the 
method so far presented. Equation (15) is therefore used to ob- 
tain the changes in loads due to specified yaw angles. A yaw of 


1° is taken as standard. Thus the change in trim will be: 


AT = Coe 


ae 





for the forebody deadrise of 22°, Figure 11(a) shows how the 
force vectors are constructed. The total load on the forebody 

is considered to be composed of two vectors each equal to one- 
half the total force acting at the center of pressure. Each 
vector acts on one side of the vee surface. For a positive yaw, 
the diagram shows that the vector on the right side has been in- 
creased by the relation of equation (15), while the force vector 
on the left has been decreased. The designation right and left 
refers to the position of vectors as viewed when looking directly 
at. Figure li and is opposite to the designation which would be 
used by the pilot sitting in the hull and looking out ahead of 
the hull. The former system is used for direct comparison with 
Figure 11. When the horizontal components are computed using the 
tangent relation of the deadrise angle, an unbalanced force to 
the left is observed. The value of 70 obtained agrees very well 


with the force of 2) that was previously shown must exist. 


= Hi, 300 1206 + 37 | z 
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The horizontal force is then 


tan 22° (21,250) = 8,570 lb. 
tan 22° (20,050) = 8,100 lb. 
Net Force 7 Lb. 


The same procedure is used for the afterbody, and the 


force diagrams on the right-hand side of Figure 11(a) show that 
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for agreement, the force due only to change in trim angle is in- 
sufficient. The value is computed as 210 lb. when from the ana- 
lysis of the static derivatives it was seen that a force of 28 
lb. must act on the afterbody. 

AC = sim Wd sear 8 
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tan 28° (4,700) = 2,200 lb. 
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Therefore it is obvious that some other force must be added to 
provide satisfactory agreement. References 11 and 12 contain pro- 
files of the transverse wave shapes for a vee-planing surface. 
Unfortunately, no empirical relation exists for the transverse 
Wave as is the case with the longitudinal wave profile. There- 
fore it was not possible to construct a transverse wave based on 
Condition 1 of Table IVY. Since no conditions in References ll 
amd 12 exactly duplicated Condition 1, Figure 27 of Reference 11 
was chosen for the closest agreement. By consulting Figure 27 of 
Reference 11, it will be seen that the point where the afterbody 
center of pressure acts is located somewhere between one and two 
beams aft of the step. It will be observed that the lateral dis- 


placement of the hull at this point results in the fact that the 


= 90% 


wake is not symmetrical with respect to it and the water level 
(not yet disturbed by the hull) is located at different heights 
above the chine on the port and starboard sides. This fact can 
be simply interpreted as the mean lateral inclination at the 
water surface. Because of this inclination, an added horizontal 
force will be created which will add to the force already ex- 
isting due to the change in the effective angle of trim of the 
afterbody. 

In the present work, a short-afterbody hull is being con- 
sidered at a small angle of yaw, and the effect of the wake has 
been found to be svabilizing. If the angle of yaw had been 
large enough for the afterbody to ride over the wave crest, the 
effect would have changed from a stabilizing to a destabilizing 
ene, likewise, if a long-afterbody hull had been used, it would 
have been found to ride on a roach, so that a small displacement 
would have moved it towards the holiow which forms on the side 
of the roach at this distance from the step. Again the effect 
would have been destabilizing, and the directional stability 
difficulties with the long-afterbody hull are well-known. 

Peele angle of yaw, the lateral displacement of the hull 
centerline at the location of the center of pressure, (i.e., 
menor fu.) att of the C.G. is 


eee 
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or ,035 of the forebody beam. At this small distance from the 


centerline, the wave slope is too small to be measured on the 





diagrams of Reference ll. Therefore the procedure shown below 


is used: 
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The angle of the wave is measured at the chine where it is found 
to be On Ge by average between the values at one and two beams 


aft of the step. The desired slope is then computed to be 


BZ peciecamedlions (237) | 


O22 


afterbody beam 5y63 


The afterbody lift force of 8,950 by acting on the slope 
of .O422 adds a lateral force of 

8,950(.0h22) = 378 lb. 

The test conditions for which the diagrams are given in 
Reference 11 do not correspond directly to the conditions of the 
hull analyzed herein. The following table gives the compari- 


son: 
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Nearest Conditions Conditions 
from Reference 11 of Present Analysis 
oneal rele 
B = 30° Semen 
CA = .73 CA = .648 
Cy =O C, = 2.63 


The number of test points in Reference 11 did not permit exact 
extrapolation, and certain bold assumptions had to be made. It 
was assumed that the differences in Cap and in trim C were 
small enough to be neglected. It was further assumed that the 
transverse slope of the wave near the centerline of the wake is 
proportional to ine angle of deadrise. This left only the ef- 
fect of the speed coefficient to be extrapolated. It appeared 
from the study of the diagrams of Reference 11 that the slopes 
are proportional to 0.9 of the ratio of the Cy involved. 
bince the lateral force is proportional to the slopes, it be- 


comes in this case 


1 we 
78 | —————— | = = 22 ‘ 
B) (a3 i : 30 ee0 


OS) 
The correct horizontal force is then the addition of 210 lb. plus 
220 lb. or a total force of 30 1b. which is now in good agree- 
ment with the force of 428 lb. which was shown had to exist in 
order that the static stability derivatives quoted should exist. 
Static Directional Stability 


It is now possible to obtain the stability derivatives 
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knowing the forces created on the hull during a i angle of yaw. 
The two forces are added together for the total transverse force, 
made nondimensional, and therefore as such represent the slope 
of transverse force curve vs. yaw angle for small angles of yaw. 
The moments are computed from the transverse forces acting at 


the arms = and X made nondimensional, and thereby repre- 


2 3 


sent the slope of the yawing moment vs. angle of yaw curve. 
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-.00291 per degree 
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-.167 per radian 

It will be observed that the above values of the static 
stability derivatives have been obtained on the basis of hull 
geometry, the empirical relations for the lift of planing sur- 
faces, and the data of the wake survey without any reference to 
the direct stability tests. It will be of interest now to com- 
pare these values with the values obtained by direct test in 


Reference 7. 





OY'/OB ON'/OB 
By computation 2 2uel -.167 


By direct test 0229 See 2 


The correlation now appears complete and for the condition 


of the hull shown the mathematical method of determining the 


static stability derivatives gives excellent agreement with test 


results. 


ir. 


A summary of the method now follows: 


For any condition specified for trim, speed coeffi- 
cient, load coefficient, and trimming moment, the 
hull is analyzed for the vertical forces acting on 
the forebody and afterbody and the arms X, and d, 
from the position of the C.G. 


Aa 


angle of yaw is given the hull and formulas (1) 
and (15) are applied to solve for the horizontal 
forces existing on the forebody and afterbody as shown 
in Figure ll. 

The horizontal force on the afterbody is corrected by 
considering transverse wave shapes from References 11 
and 12 and the additional force these wave shapes add. 
With the horizontal forces actingat X, and h , the 
static derivatives may be computed for the le yaw 
angle, giving the slopes of the desired curves per de- 


gree yaw. These may then be converted to radian mea- 


Sure. 


The mathematical treatment is not complete until the ana- 
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lysis has been made which produces the expressions and equations 
necessary to convert the static stability derivatives into dy- 
namic stability derivatives. 

Rotary Stability Derivatives 


Returning to equation (2), the moment N may be expressed 


; oY +N ae 
™ na | Ae 3 hz 7 


which in reality is an identity for the yawing moment N and 
yields 

ie th) ee. 

However, equation (16) is used to derive an expression 
for the rotary stability derivatives by giving the hull in Fig- 
ure 1(a) an angular velocity "r" about the Z-axis. If the angu- 
lar velocity "r" is referred to the C.G., such an angular veloc- 
ity will create a linear velocity ENG at the forebody center 
of pressure and a velocity ri, at the afterbody center of 
pressure. 


By def inataon: 


N= NioB eyo (17) 
end 
Y= ¥',8 ByMp* (18) 


The effect of giving the hull an angular velocity is to 
add or subtract from the forces resulting from an initial yawing 
angle another yaw angle AB which for small angles on the forebody 


has the value wey and on the afterbody rhifV . The effec- 
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tive angles of yaw are then: 
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Substituting equations (17), (18), and (19) into equa- 
tion (16) gives 
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where AN is defined as 
a 1) oy2 (: Oe 
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Using dimensionless coefficients: 


i= ! ! 
N N g BY + AN 


The resulting moment on the hull is the sum of the moment 
due to initial yaw angle plus some change due to the angular ve- 
locity. For the purposes herein, interest is concentrated in 


the AN' created by the angular velocity, that is, for zero ini- 





tial yaw angle. 
. AN | “Kata, (a) ~ Aang (EP) 
Neh, 
sng le) l8 
Ne 


Differentiate with respect to the dimensionless parameter r' = 
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This represents the equation necessary to determine ana- 
lytically the rotary derivatives, knowing the static derivatives. 
The quantities Ne and Nis are determined for the planing con- 
ditions, assuming that for small disturbances from a straight 
course the forces FS and Fy, continue to act at the centers 
of pressure of the forebody and afterbody. 

Using equation (1), 

ao ee 


and also 
N+ iy \yy-nN 
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This again is an identity in Y 3; however, if once again 


equations (17), (18), and (19) are restorted to, the following 
is obtained: 
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Making results nondimensional 
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_ The same situation exists for transverse forces as for 
moments. That is, the total Y force is the sum of effects 
from initial yaw angle and added angular velocity. Interest 
is again centered only in that transverse force due to angular 
velocity. Therefore taking the expression for AY' and differ- 


entiating with respect to the dimensionless parameter r= 2 


r ag Lo 


Y' =-N', . (2p) 
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ANALYSIS FOR CONDITIONS 2 AND 3 


The analysis of the hull for a speed coefficient of 2.62 
was a straightforward application of the equations for lift of a 
prismatic vee-planing surface, the equations for center of pres- 
sure measurements found in Reference 10, and the shape of the 
longitudinal wave profile as computed in Reference ll. OQne of 
the limiting conditions in Reference ll is that the speed coeffi- 
cient be no less than 3. Difficulty was encountered when the 
hull was analyzed for speed coefficients of 2.46 and 2.3. The 
value in equation (10) of H, contributes the largest portion of 
the computed wave height. The following values of the trim and 


speed coefficient as entering variables for the three conditions 


are noted: 
a e all mote <o el 
Condition C at ae C C/G 
V V V 
i 12.6° 2.62 5.88 2.366 1.783 Bee 
2 1 2a olka ILarAle 3.16 
3 10.5° 2.3) 5.18 2.225 1.665 Bish 
It has been mentioned previously that at a ©, of 2.62, 


V 
the formulas for the longitudinal wave profile were used beyond 
the range of the wake tests which went down only as low as 

Cy, a. oO 5. at Cy = 2.16 and 2.34, this extrapolation had to 
extend so much further that the trial applications indicated 


that it was not possible to obtain the proper wetted length of 


the afterbody by this means. In absence of the actual wake data, 
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therefore, the following alternate procedure was adopted. It was 
clear from the basic theory that the forebody load will be de- 
creasing and the afterbody load will be increasing with a de- 
crease of the speed coefficient Cy - In Condition l, the fore- 
body load was found to be 81% of the total load. It is assumed 
now that the forebody load becomes 75% of the total load for 
Condition 2 and 70% for Condition 3. Using the trim angle data 
obtained from the test and the planing relations formula, the 
forebody wetted length and the position of the center of pres- 
sure were computed. The afterbody load was taken as the dif- 
ference between the total and the assumed forebody load. The 
necessity to obtain the balance of forces and moments therefore 
dictated a definite value for the distance from the C.G. tne 
position of the afterbody center of pressure. A suitable set of 
values for the afterbody wetted length and afterbody effective 
angle of trim was obtained by trial and error, so that the fric- 
tional resistance of the afterbody was properly accounted for in 
the balance of moments. In absence of the available wake data 


for the low values of C the above procedure can be considered 


V 3 
as the best substitute. 
It was assumed that the transverse wave profiles can be 


extrapolated to the low values of OC. in the same way as it has 


V 
been shown in detail for Condition 1. With the mean lateral 
slope of the wave obtained in this way and with the values of 


the afterbody angle of trim as obtained as explained above, the 
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changes of the effective angle of trim on the two sides of the 
hull and the resultant lateral force were readily computed. The 
results of these computations for Conditions 2 and 3 are given 

in Appendix II. A recapitulation of all the hydrodynamic planing 
characteristics is contained in Table V for the three conditions 
considered. Figure 9 contains the afterbody wetted areas for 
each of the three speed coefficients investigated for comparison 
purposes. Figure 11 contains the hull force diagrams for the 


three conditions investigated. 
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MATERIAL AND TEST PROCEDURE 


The model used in the tests was that of a conventional 
large flying boat designated as Stevens Model 06. The model di- 
mensions appear in Tables I and II. A diagrammatic sketch in 
Figure 2 shows the essential parts of the standard seaplane yaw- 
ing apparatus which was prepared using Reference 13 as a guide. 
The carriage on which the model was mounted allowed for freedom 
in trim, heave, and yaw. Freedom in heave was necessary to allow 
for correct wetted length although the heave for each run was 
not Peeeried. The apparatus was equipped with an unloader which 
permitted any desired load on the water to be obtained regardless 
of the weight of the model. The yawing tests in the present 
work were run at a constant load coefficient and at zero heel 
angle although it was possible to vary the heel angle. Angle 
of trim and yawing angle were recorded for each run. 

No apparatus existed in Tank No. 1 for the direct deter- 
mination of the sidewise force due to yaw angle. Therefore an 
indirect method was used which proved fairly satisfactory. The 
moments about the C.G. and about an arbitrary point 3.56 in. 
forward of the C.G. were computed at the same fixed trim which 
coincided with the free~to-trim angle obtained for zero yaw 
angle. The moment about the point forward of ©.G. minus the 
moment about the C.G. is equal to the sidewise force acting at 
the C.G. times the distance between the two points. This equa- 


tion may be made nondimensional and differentiated with respect 


ee 





to B , which will lead to the following result: 


N., = Ny =- ¥ xX 
N,' ~ N,' = an 
ON,' ON,! 
OB O8 dy 
—— Se 
where 
N = dimensional force about the C.G. 
N,! = nondimensional force about the C.G. 
N, = dimensional force about forward pivot 
Ne! = nondimensional force about front pivot 
X = dimensional distance between pivots 
X' = nondimensional distance between pivots 


Therefore, the model was equipped with two sets of pivots, the 
distance between the pivots being 3.56 in. The model was bal- 
anced around the rear pivots where the C.G. was located 2.00 in. 
forward of the step and 6 in. above the keel. This corresponded 
fora tui scale C.G. position of 3.67 ft. forward of the step 
and 11.0 ft. above the keel. The position of the model inthe 
Carriage is shown in Figure 3. 

It was necessary to make three runs to record arene 
data. The first run in the free-to-trim attitude made use of 
the rear pivots of the model where the C.G. of the model was 
established. A standing yaw was recorded, the model sent down 


the tank at a predetermined speed, and the running yaw recorded 
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as well as the trim. The difference in the standing and running 
yaw was multiplied by the spring calibration to give the moment 
in inch-pounds (calibrated in inch-pounds per degree change in 
yaw angle). The resulting moment was made dimensionless and 
plotted as yawing moment vsrunmng yaw angle. The springs which 
restrained the yawing motion of the model were available in four 
varying degrees of stiffness. The calibrations for the two used 
are shown in Figure 12. These two springs were distinguished 
by red and yellow paint. A dashpot was also’included which 
damped out oscillations in yaw. The spring was located 15.5 in. 
from the yawing pivot. 

Reference 2 contains a good description of the procedure 
sometimes followed when hulls exhibit directional instability. 
The method consists of actually pushing the model while running 
in order to get points on the curve which would otherwise be 
impossible to obtain due to a combination of instability and 
spring stiffness. Figure shows the relation between standing 
and running yaw. Since the yaw angle could be read only to the 
nearest 0.1, springs were interchanged at various times to im- 
prove the scatter of the points on the curve. 

Once the running trim angle was determined, the model was 
locked in trim at the free~to-trim angle for zero yaw and run 
again for angles of yaw from 6° port to 6° starboard. Immediate- 
ly following this, the model was locked at the same angle of trim 


on the front pivot and runs were again made covering the same 
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range of yaw angles. The trim lock was of simple friction type. 

All the testing was conducted in Tank No. 1 of the Ex- 
perimental Towing Tank, Stevens Institute of Technology. The 
tank is 95 ft. long and of semi-circular cross section with an 
added dock section 6 ft. long to facilitate preparation and 
mounting of the model. The carriage is suspended from a mono- 
rail which follows the centerline of the tank. Power is re- 
ceived from a towing line which is driven by a mulitiple-diameter 
cone pulley. Provisions are made for further varying speed by 
interchangeable gear boxes. Power to the towing line is supplied 
by direct current until the model is up to speed at which time 
an A-C. synchronous motor cuts in to operate the carriage at a 
constant speed. 

The water level of the tank had to be carefully control- 
led so that the carriage would have no restriction in heave. 
The water level for fixed trim was generally lower than for the 
ree-to-trim case. The criterionof correct freedom in heave was 
that the carriage horizontal balancing arm should extend slight= 
ly above the horizontal position when the model was at rest in 
the deck. 

Prior to running, ali three test conditions shown in 
Table IV were checked with Reference Ui in order to ensure that 


the hull was not in the porpoising range. 





‘TEST RESULTS 


Tables VI through XIV tabulate the results of the tank 
yawing tests, and these results are shown as moment vs. yaw 
angle curves in Figures 13 through 21. Figure 12 shows the 
spring calibration for the two springs used for determining the 
value of yawing moment. Table IV gives the three conditions 
used in the present investigation. Figure 10 shows a photograph 
of a typical run as obtained during the tank testing program. 
Table XV and Figure 22 are considered to constitute the final 
results, for, in these, a comparison is made of both the computed 
and test values obtained. 

The three test conditions chosen correspond to conditions 
which were predicted to exist just beyond the unstable range. 
That is for a Cy less than 2.31, the model becomes unstable as 
shown on a composite static stability chart for a similar typical 
hull configuration. The test results show that with a decreas- 
ing speed coefficient, the hull becomes increasingly stable. 
However, the stability is represented by three stable branches 
of the curve with discontinuities existing in the vicinity of 5° 
port and starboard. For the lowest speed coefficient, the hull 
was the most stable, with the discontinuity between the three 
stable branches existing at about 3 port and starboard. It was 
predicted in the mathematical analysis that a decrease in the 
speed coefficient would lead to an increase in stability. How- 


ever, the stability predicted for the lowest speed coefficient 
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was approximately one third less than the static stability found 
by test, in the case of yawing moment vs. yaw angle. 

All the free-to-trim moment curves for the three conditions 
check favorably with similar curves previously reported (Refer- 
ence 2). Difficulty was encountered in setting the model at 
zero yaw angle. No accurate apparatus exists in Tank No. l for 
such a procedure. So much water had to be let out of the tank 
that the measuring device used had to butt against the curved 
sides of the dock. The red and yellow springs were constantly 
being changed in order to obtain better curves. Each time the 
spring was changed, the point of zero yaw was shifted. The fact 
that the zero yaw calibration was not good has shown up in the 
moment curves as a horizontally displaced curve. The apparatus 
was set at zero heel but it was noted that quite a bit of play 
existed in the pivot support which was equivalent to inducing 
an angle of heel on the model. Figure 18 is the only curve 
which shows a vertically displaced axis due to some induced heel 
angle. The displacements of the axis in either direction for 
these tests were considered small enough to have no effect on 
the slope of the moment curve for Zero yaw angle. 

Concern in moment curves for the fixed trim case was cen- 
tered in small angles of yaw. Since the change in trim angle was 
not too large for small angles of yaw, fixed trim was taken at 
the value of trim for zero yaw angle. Using the method already 


suggested of subtracting the slopes of the moment curves for 
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front and rear pivots at the same fixed trim, the following are 


the solutions for’transverse force vs. yaw for the three condi- 


tions tested: 














oy = 2.63 
ON," 
Front pivots = = ,0075 
OB 
N ! 
Back pivots — = = ,0050 
p 
oo yor, 
eo B © B oe 
_ OY! _ - .0075 + .0050 
OB fobs 
ou. ,00383/degree 
= ,219/radian 
Cy Z110 
)N,' 
Front pivots ra = -,0052 
Q Ny! 
Back pivots SEP = =-,0020 
_ OY! _ = .0052 + .0020 
® B 2053 
ot .00L9/degree 
= ,281/radian 
Cc 
_ oO N,! 
Front pivots 5 : = -.018 
dN 
Back pivots <== = =. 005 


Qo PB 


A 9 a 
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7 Oy! _ ~ .018 + .023 
o8 053.” 


St .00766/degree 


2 


= ,439/radian 








MISCELLANEOUS 


The testing program in Tank No. 1 began with Condition l, 
wherein the speed coefficient was 2.63. This was oe case where 
the best agreement was desired in order to form a basis for the 
subsequent analysis. Unfortunately, due to the lack of previous 
experience with the apparatus and tank testing techniques, the 
experimental data in this first condition were obtained witha 
larger amount of difficulty and with less precision than in the 


subsequent work. Consequently, Figures 13, 14, and 15, which 


correspond to Condition 1, show quite a scatter of points. With ’ 


the existence of such a scatter it was quite difficult to fair a 
curve through the points. The general trend of the curve seems 
obvious but the reported vaiue of the slope might have been 
better. It is significat that similar runs for the lower speed 
coefficients contained much less scatter, which may be partly 
due to the true effect of speed and partly due to the increased 
famitiarity with operating techniques as the testing progressed. 
t 1s regrettable that sufficient time was not available in which 
to re=-test Condition 1 in the hope of better agreement with pre- 
viously reported static stability derivatives for this condition. 
During the entire testing program the yaw scale could he 
read only to the nearest 0.1 of the yaw angle. For Condition 1, 
a change in 0.1 of the yaw angle meant a change in the dimension- 
less yawing moment of .005, which, for the scale used, made a 


Certain Sameuny of seaupter almost inevitable. 
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As pointed out previously, the continued horizontal dis- 
placement of the axis was due to difficulty in resetting zero yaw 
when changing springs. The vertically displaced axis which oc- 
curred in only one case was in all probability due to excessive 
play in the pivot support. 

The validity of considering slopes for front and rear 
pivots at the fixed trim indicated for free-to-trim at zero yaw 
seems adequately illustrated by the agreement between computed 
and calculated values. Therefore, since only small angles of 
yaw were considered, the differences in the slopes for fixed trim 
at zero yaw were satisfactory measures of the slope of the trans- 
verse force at zero yaw. The slope of the transverse force vs. 
yaw angle was considered linear for the range of small yaw 
angles. The only satisfactory measure of the moment curve vs. 
yaw angle slope was obtained at zero yaw angle for the free-to- 
trim condition. The moment curve is evidently much more sensi- 
tive to trim angle than the transverse force curve. 

In the early stages of the testing schedule it was felt 
that water interference with the tail cone might possibly be con- 
tributing to the scatter of points as observed in Figure 13. In 
order to.check this effect, a check test with spray or breaker 
strips was made. Subsequent to the installation of these, data 
were discovered in a previous test for tail cone spray strips on 
a Similar model. For the speed ranges considered, it was found 


that previous results did not contribute anything to the data 
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and since nothing was known conerning the best position for the 
tail strips, they were removed for subsequent testing at the two 
lower speed coefficients. The addition of the strips, it was 
felt, might have produced more uniform results. 

The operation of the apparatus depends on the fact that 
the model and retaining springs form a stable system. Pierson 
(Reference 2) has suggested a method for trying to obtain points 
on the moment curve which normally would be impossible to get 
due to the stiffness of the spring. Figure shows a typical 
yawing moment curve with the equilibrium position that the model 
will assume depending on the angle of standing yaw. The process 
involves pushing the model while in motion to force it to assume 
an equilibrium position other than what it would normally take. 
The process has been used to advantage on unstable hulls. In 
the case of the lower speed coefficients where the hull has three 
stable branches, this method was tried but without success. 

The type of mechanism used to hold the hull at fixed trim 
was probably not as good as it might have been. The friction 
lock was too close to the pivot point of suspension. It was 
originally designed for lock on the front pivots only, until it 
was decided that the moment curves were best compared at the ex- 
act same angle of trim and therefore must lock on the back pivots 
too. The friction lock was attached to a bracket on the midship 
section and to a projection on the pivot support rod. When the 


model was set at an angle of trim, it is believed that because 
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the apparatus was not correctly aligned, a slight angle of heel 
and angle of yaw were set into the model. This can account for 
the difficulty of zero yaw and zero heel for fixed trims. 

The choice of the model used during the testing position 
of this investigation was dictated from consideration of the 
wealth of previous test data already in existence for this par- 
ticular model. Dueto personal inexperience with model testing 
technique, it was felt that use of such a model would provide 
many desirable check points. Unfortunately, the existence of so 
much available test data also meant that the model hed received 
an excessive amount of use and could not be put into very good 
shape. It would, in all probability, have been much better had 


a newer model been immediately available. 
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CONCLUDING REMARKS | 


The methods described in this report permit the lateral 
forces acting on a forebody and afterbody to be estimated by 
equations (14) and (15) on the basis of the yawing moment and 
lateral force obtained in a yaw test in a straight towing tank. 
The rotary derivatives are then obtained from equations (20) and 
cal). 

Alternatively, it is possible to base the calculations 
on the angle of trim obtained from a simple zero yaw test, and 
to make an estimate of the wetted areas and centers of pressures 
on the basis of the empirical formulae of Reference 10. 

Figure 22 shows a graphical comparison of computed and 
test values of the static stability derivatives. It will be ob- 
served that the greatest divergence between tested and computed 
values of QN'/QB occurs at the lowest value of speed coeffi- 
cient used. The reason for this is not too clear. Down to 
C. = 2.31, the hull has been rapidly becoming more stable with 


V 


decreasing C From the composite static stability chart of 


V e 
Reference 2, it may be seen that with a further decrease in Cy 


from a value of 2.31, the hull becomes unstable. This indicates 


that values of C, from 2.3 to 2.4 lie in a transition region. 


y 


It is possible, therefore, that there is some change in the trans- 


verse slope of the wave in this region which affects the afterbody. 


Such a change in the wave form would not be predicted by the ex- 


trapolation used in computing the lateral force due to the trans- 
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verse wave. 

The procedure developed in the analysis section of this 
report for predicting the static stability derivatives shows a 
reasonable agreement for the speed range considered for hulls 
of conventional design with results obtained from standard yawing 
tests. Unfortunately, it was not possible to check agreement 
in regard to the rotary derivatives since no experimental infor- 
mation exists concerning these. 

It 1s believed that a more complete comparison between 
tested and computed values of the derivatives was not possible 
and a closer agreement was not obtained because of the following 
conditions: 

a. Absence of data on the longitudinal wave profile at 

low values of Cy : 

b. Insufficiently complete and accurate data on the trans- 

verse wave profiles. 

c. Absence of sufficient experimental data on the rotary 

stability derivatives. 

It is suggested that steps be taken to correct the above 
deficiencies, and that the same type of investigation as described 
herein might profitably be applied to other types of hulls and 


at a wider range of speed coefficients. 
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APPENDIX I 


It is desired to obtain the equations of motion for a 
body moving in a nonviscous and incompressible fluid. From any 
standard text on vector analysis, one finds the following ex- 


pression for the linear momentum and the angular momentum: 


no (Y) = y (Sv + Ox wy 
dt : dt e 
fixed moving 
system system 
dH = dH + OxH 
dt ‘ dt : 
fixed moving 
system system 


where the vectors have the value 
= uli + vj + wk 


pi + qj + rk 


my DO) <4 


ho eee hk 
x Vy Z 
Evaluating the vector product by the matrix, equating like compo- 


nents, and applying Newton's law to the coordinate axes fixed in 


space gives 


1 J k 
—> —p 
Lx MV= |p q r 

! ! 

ue mV Mm. W 
ey a m,'u = m,, Vr as ro = X 
E, = mM, !V + m, 'ur - Ue le a 
Ps ~ ey f mM, ' VP ~ m,‘uq = Z 
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M (at =F (resultant force) 
dt : 
fixed 
system 
~F = Poi + Fed + Pk 
1 3 k 
> > 
0) ise q r 
h h h 
x Vy Z 


i(h,g ~ rh.) + j(rh, - ph, ) 
° k (phy eae) 


ee 
(gg = G (resultant torque) 
fixed 


system 


Sac ° e 
G = Gj = God os oe 


G, = nh = oe 7 qh, = L 
G, = By + rh, ~ ph, = M 
G, = h, + ie = qh, = . 


The forces F are the externally applied forces acting in a 
Newtonian frame of reference, the components of which are found 
by equating like components of the original equations. From 
geometrical considerations, one considers an element of mass 
“dm" and determines the contribution of this mass to the angular 
momentum about the three coordinate axes. It is shown that 

h_ = Ap - Fq - Er 


Xx 


h 
ME 


-Fp + Bq - Dr 


= 60 = 





h, =s—eip = Dg + Cr 
where 
hea] D= I 
x Vz 
B= I B= I 
y XZ 
C=] F= 1 
Z, XY 


which for a symmetrical body taken along the principal axis, the 


products of inertia will vanish, leaving 


hy = Ap 

h = Bg - Dr 
a 

h, = - Dq + Or 


The problem to be investigated is concerned only with 
yaw; therefore, it is supposed that the motion in the horizontal 
plane may be considered separately. This leads to the more 
simple form of the equations of motion since roll and pitch and 


heave are made zero. The final equations desired, then, are as 


follows: 
m,'u ~ m, ! Vr = X 
my!V +m, ‘ur = Y (I-1) 
r= 


REF. 
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In order to solve these equations, one makes several 
transformations by putting the various terms into dimensionless 
parameters. With reference to the accompanying figure, the 
following substitutions are at and put into the original equa- 
tions: 

wu = ¥ecos 8 

v= - ¥ sin B 


lies p — ) Gam p B 


<de 
{ 
) 


V San pt ¥ B cos B 


_ Ge 
rt” at 8 
Se Wa 
= t . ! , = 
m,'V cos B mV B sin p + Mm, NO serile oo aaen 
-m,! V sin B - mV 8B cos Bp + m,' V @ cos B = Y (I-2) 
I. 6=N 


When the original equations were developed, the mass was 
merely referred to in three different directions coinciding with 
the coordinate axes. iceman the aid of the "ellipsoid approx- 
imation," the mass in the X and Y direction have the fol- 


lowing standard seaplane notation: 


3 
fl 


m + mk 
m,! = nie mK5 
where m is the mass of the seaplane, Mm is the mass of the 

displaced fluid, and ky is the virtual mass coefficient in the 


X direction. The following dimensionless substitutions are 
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then made: 











X! = X m 1 
Pyeye it Py 
E 
m + Mek 
— fe 5 2 mee 0 > : 
ae 2 
N I, + K'Ig 6 
q' 7 fc 3 aaa Be 
a” 2 
m,bV cos B m,VbB sin Bp m,Vb@ sin B 
eee en See a” | y, 
ve ve ve | 


A quantity s 


_ Vt 
ee 


is introduced which is a distance defined as: 


and therefore all derivatives may be replaced in terms of this 


dimensionless parameter as follows: 





da _d ds 
ait ds@dt 
aiVa 
dt b ds 
dC). ( yo = BAC) 
ds V dt 
b 1b db 
gi-— , ete 5 gt = =r 
,.d'iv,ri a 
at % peas 
ane oe 
be be 


Substituting these into equations (I-2) and limiting the 


discussion to small values of 8 , it follows that the three equa- 


ioe 





tions take the following form: 





mV! 
=  - mVB'p + mor 'B 

- mV! 
Y! = —~— 8 - mB! + myr' 
N! = nr" + NE Da 


For steady course, it is noted that V!' is at least first order 


in 8 and r' since X! is zero for steady course, and further 


that V' is at least second order in Y' and N'! so that V'! may 


be made equal to zero. Since the motion under consideration is 
concerned only with the horizontal plane, the force Y!' and the 
moment N' may be considered to be functions of B and r' , 


Expanding into a Taylor series and neglecting higher order terms 





gives 
= ner + i'r 
oY! 9 Y! 





il 
Ke 
“A 
+ 
" td 
“TH 
“TD 


For initial values taken as zero, changes are small and assumed 
to remain in the linear region. 
By substituting these expressions, the following differen- 
tial equations are obtained: 
I ! - ! s ! = 
m8' + X gf r (m, vg ane 0 


nr" - Nt ort -wyN = 0 


r! ge 
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Using operator notation, two simultaneous equations are 
solved; the solution is obtained from the quadratic formula. The 
solution represents the roots of the quadratic and is a solution 


to the differential equations. 


= 


@| 2 


B(m,D as ¥'y) ~ r'(m) ae Ee a0) 


Ei ia la) = N',B = 0 


m, - Y'L, ; (nD - Rea 
m,D + Y' N! 
e B p 
aC ~ 
m nD” + D(n¥' - mN'.,) + ee - Nt ,(m, = voy =41G 
D- 01> 
; ae f = mN',) 
em, 


e 
2 
{or ee [r te a rh) | 


emon 
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APPENDIX II 


Mathematical Calculations for Condition 2 of Table IV 

For Condition 2 it was found that the speed coefficient 
was much too low for laying out the longitudinal wave profile 
as shown in Reference 8. The process then was one of trial and 
error in order to determine forces and distances as shown in the 
mathematical analysis section of this investigation. For the 
lower speed coefficients it was impossible to get the added in- 
crease in afterbody load that was necessary since the load car- 
ried by the forebody decreased as the speed coefficient decreased 
and therefore the afterbody load must be increasing. 


Forebody Conditions 


C= Ialeo 
C= 2.6 
2 = 
(C.,) 6.0) 
CA = 0.6 
BAGS). 
CT ema te? 
B 
Cr = ,262 - .143(.448) = .262 - .06) 
B 
=alog 
C, = .262 
O 
om 
ari = ,01832 


a 
= .0120(1.51)2 + .001572(1.52)* = .01476 + .00357 


= .01832 


= 66% 





ees od 


> 
l 


£6062,3) (1000) = 3'7,),00 


_ 37,400 
(Ap = Gos 11.20 * 38100 


ri 
C = 75 ae = ~669 
P Con 
3.06 Ca) + 2,h2 


P = ,669(15.1) = 10.11 
A, = 10.11 - 3.67 = 64h 
Afterbody 
(A), = 49,750 - 37,400 = 12,350 
For b = 6.8, 
po uum coculigs 0 = 13 


For 


Lg ~ 1,937(1935) (46) 





C, = 216 - » 182(. 399) 
B 
= 13h 
C.. = ,216 
O 
oy 
re 01687 . 
= ,012(1.47)2 + .001068(1.47)* = .01457 + .00230 
= ,01687 
eS. Ley 


e Z 7 
C, = .75 - IGES +4 Im = .75 - ,056 = .69) 


ee 1.47(6.8) = 10 ft. = 120 in. 
f 


a 2 20 in. 


= 67 = 





a ae 


Pe="Foord0) = 6y9N ft. = 83.2 an. 


= 1 = 
tan >} ee 532 





nm .1796 
3} = 28° 
Ly = 13.83 + 3.67 (Figure 9) 
Sel 5 Omer 


Forebody Frictional Resistance 


_ 943 (SS 


S 
(Re) = 5.99 x 10 
Mm 1,08 x 107 
_ eer 
(CE) = pI? 200516 
_ i Z 
(Ca) -00516(1.51) (359) = ,008)2 


R = .008)2 ‘ae (1935) (100) 


1578 lb. 


Afterbody Frictional Resistance 
9.43 (2922) 
eee - — = 3-97 x 107° 
m 1.08 x 10 





— ,OO7 _ 
(Ce), = 77 * .00562 
(Ce), = ,00562(1.1,7) (aa = ,00937 


R = .00937 (10937 | (1935) (46) 
806 lb. 


Total R = (1578 + 806)1.25 = 2980 lb. 





=960-— 


Moments about C.G. 
36, 1OGtee ie = 1OCzNFe0) + 117.50(12,350) 
245,500 = 29,800 + 216,500 


a 


= 26,300 


Forebody Transverse Force (consult Figure 11(b)) 
AT =.37L (for 1° yaw) 


z ees) 2 
=~ les 2 Pa Pov a 
(4), = 19,050 (22 5 peer ls 


tan 22°(19,678) = 7.950 


Il - 


tan 22°(18,321) = 7,399 


Dole 


tl 


Net Force 


Afterbody Transverse Force 








(A), = 6,175 ccna = 6.159 
(A), = 6,175 (ew) = 5,890 
tan 28° (6,459) = 3,36 


tan 28° (5,890) = 3,13h 





I 
Lo 
© 
NO 
es 
Oo 


Net Force 


Additional Afterbody Transverse Force (Figure 32 (Refer- 


ence 7)) 


distance from ¢ =e 5 
slope 1 beam = 9” 


slope 2 beams = 16.3 
Average = 17.9 








= FO 4 


tan 17.9° = .323 
ae 50 


center of pressure travel = =—— = .305 


5763 
vertical distance = (.305)(.323) = .0985 


WIS 


total = 2(.0985) = .1972 


wave slope = ele = 3029 


added force = (.029)(12,350) = 3.58 lb. 





Correction: 
ili Oe _ 
356 sont ; 55 = 290 lb. 
Zain 
OY! _ 2(551 + 592) 


B 1.937(1935) (100) 


.0061/degree 


@ 


ul 


»349/radian 


ey 


Nt _ 21551(6.4k) - 592(17.50) | 
1.937(1935) (1000) 


| 


Q 
mee 


- ,0036/degree 


= - ,208/radian 


Mathematical Calculations for Condition 3 of Table IV 


Forebody Conditions: 


C= “Toes? 
ye ores 
(G5) ¢ =" 507 
CA = 0.56 





Cem! = 29 
8 54h 7 


rr? 
2 





C. = .270 ~ .143(.456) 





B 
= ,205 
are 
0 
- 
wT = 90203 
=nOmeO We29ey) + .001737( 2.78) 
= 01551 + .00483 
= 0203) 
ah = or 
A= .56(62.3)(1000) = 34,900 





“(eye ewe 3 ob 208 
cos 10.5 
Can eee 
Ps a7) + 22 
= ,052 
P = ,652(16.7) = 10.9 


> 10.9 = 3.67 = 7.23 


Afterbody Conditions 


(A), = 49,750 ~ 34,900 = 14,850 
For b = 7.5 
ee. «1595 


B 1.937(1720) (56) 


Cy aa 0235 ae » LOZ aa > 
p 


= 0235 ia: -Omo 


ae ey 


-~ 7) 








ft 
O 
For Tt = 10,3 
ie 
et =OnnO> 
= ,0120(1.246) + .001296(2.h) 
= ,O19) + .00311 
= 01805 
h = 1.55 
porate (ied) 
= 11.62 ft. = 139.6 in. 
AY = 2. = 279.2 in. 


= = ij — = a 
Cy = on = eG? Feoeio 75 B07 1 0679 


P = .679(11.62) = 7.9 ft. = 94.7 in. 


A 5 = Ulw9 + 3.67 = 15426 


Forebody Frictional Resistance 


6.908 (263"| 





(Re) = ~ = One x 10° 
a0) se ile 
_ 007 _ 
(Cy), = Domne = 00522 
=| al = 
(Ch), g00512(1.67) 537 200922 
R = 00922 (2937) (1720) (100) 


eeeb. 


I 





= 73 = 


Afterbody Frictional Resistance 





8.988 1162| } 
(Re) = = = lh.40 x 10 
i | 1 Oeee 10 
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Fim 1.35 
(Cp), = -0055(1.55) -i:5] 
= 00965 
R= 00965 (2.932) (1720) (56) 
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Mews Re= 2,25(1535 #7900) = 3,050 


Moments About C.G. 
35,500(7-23) = 14,850(15.26) + 10(3,050) 
257,000 = 226,000 + 30,500 


#) 
= 256,500 
Forebody Transverse Force 


17,750 20.5 + 2304) = 18,382 
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R 10.5 
(A), = 17,750 peal = 17,108 
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tan 22 (18,382) = 7,426 


tan 22°(17,108) = 6,912 
Net Force = 51) lb. 


Afterbody Transverse Force 


(4), = 7,25 fous sat =a Tol 
10% 3 
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(A), = 7,425 (20.3 = <tr) = 7,085 
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tan 28 (7,76) = 4,130 


tan 28 (7,085)= 3,769 
Net Force = 361 1b. 

distance from ¢ = ok 

slope 1 beam = 2h° 


slope 2 beams = 20 


I 
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NO 
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Average 
tan 22° = .hO 
center of pressure of travel = see = .266 
vertical distance = (.lj0))(.266) = .1076 
total = 2(.1076) = .2152 
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wave slope = = .0287 
Te5 
added force = (0287) 14,850 = 26 
Correction: 
1 Coe. 
125( sh —) ao 328 
Pas 


Oy! _ _ 2(514 + 689) 


OB — 4,937(1720) (100) 

= ,00722/degree 

= ,)i)/radian 
ON! _ 2[514(7.23) - 689(15.26) ] 
QP 1.937(1720) (1000) 

= -,00)08/degree 


I 


-,23)/radian 





TABLE I 


MODEL PARTICULARS 


Stevens Model No. O06 
Scale 22 
Beam at main step, in. 5.45 
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Angle between forebody keel and baseline, deg. O 
Angle between afterbody keel and baseline, deg. fee's 
Deadrise at keel and main step (average) Za 
Deadrise on afterbody (average) 2u 
Height of main step at keel, in. re 
C.G. forward of main step, in. eae 
C.G. above baseline, in. Gace 


TABLE IT 


AFTERBODY DIMENSIONS 
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Station, Baseline Baseline yee 
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LONGITUDINAL WAVE PROFILE 


(X and H are Fractions of Maximum Beam) 
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Conditions: GCG. = 2.62; C, = 0.8; = 12.6° 
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TABLE VIII 


FIXED TRIM; BACK PIVOTS 
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TABLE IX 
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) FIXED TRIM; BACK PIVOTS 
lie A 


i) = =O 
@ = 2.16 {= 236 
ch} = 0.8 
Ba A = 4.66(1b.) 


No Spray Strips 


i 
2 
3 
\ 
5 
6 
7 
8 
7 


See EO 
Creo ce = Os Wn WwW mM HE O 








TABLE XI 


FIXED TRIM; FRONT PIVOTS 


on ed A, = 4.36 
onc +o = 38 

eae 

va oe A = .66(1b.) 


o Spray Strips 


-.02)6 
-~.0328 
+.0309 
+0306 
+.0225 
+,016) 
+,0103 
+,00205 
+,00)1 
ant Oo ime 
7Oaie 
P0504 
+.0389 
Ou 
Oar 2 
-.0082 
-.0123 
- .00205 


1 
2 
3 
4 
5 
6 
7 
8 
2 


MP PrP YP PPP PP 
Gro GG -—l OO. Vt & Ww © F’ © 








TABLE XII 


FREE~-TO-TRIM 
C., = 2. 3h A = 4.16 
CAce 0.8 ae 78 
No Spray Strips = },66(1b. ) 


1 
2 
3 
4 
5 
6 
7 
3 
Z 


bo 
HH oO 





JY 
NO 





TABLE XIITL 


FIXED TRIM; BACK PIVOTS 


B= 10.5 A = h.36 
Gee 2,3); OO) 
Ch = 008 A = 1). 66Gb. ) 
No Spray Strips ° . 





Run | (Still), | (Moving), ay |. Ne N, a alla 
NO. deg. Gece in.-lb6. |ft.-1b. deg. 
0 O O O O O O 





Oo @OnN A WM FF W DY FF 





} 
© 















Cee Se 


hoo Oo et i 
mM fF O OO © ~w 





NO 
Lo 


i. Spring change beginning with Run 


be 
sl Or 





6.1P 
Bou 
5.0P 
SP 
4.OP 
Babe 
3.0P 
2.1P 
2 ile 
ele 
ones 
eS 
4.98 
3.98 
2.08 
1.0P 

SE 

oS 
eS 
ines 
2.08 
aes 


Gee 
o20F 
5.6P 
5.1P 
4. 8P 
Ieee 
3.9P 
Pee 
1.8P 


6S 
mie 
. 8S 
. 8S 
258 
ao 
2 3P 
oe 
ae 


ae Os 


4.68 
elie 


° se ee e ee ee ee e@ ee @#  e@ ee ee ae ee @ @ ee e@ ea e e 


ae 








Nese 

2.3 

Olas, 

o Spray Strips 


i 
2 
3 
Ly 
5 
6 
7 
8 
7 


TABLE XIV 


FIXED TRIM; FRONT PIVOTS 


oo 


4. 36 
30 


= .66(1b.) 








+8 +Q 











om +V by 
a _ _sTep | 
| 
C.R 
(a) (b) 
STATIC art - A 
SEAPLAN 


AT A GIVEN ANGLE “OF YAW 
FIG. 1 


+ 





_ DIAGRAMMATIC SKETCH OF 
SEAPLANE YAW APPARATUS 


MAIN TOWING GATE 





YAW ANGLE ADJUSTMENT 


DASHPOT 
Wy 
aa 
l| 
iw 


YAWING SPRING 


HEEL ADJUSTMENT 


FIG.2 





eos 
INSWSONVYENV SNLVEVddY SLAIdWOO 








= 


rt 

=z 

Ltd 

es 

© 

=. 

wD 

z 

z 

> 

/ SPRING LINE 
\NITIAL SETTING f ODEL 
PORT / STBD 

| , | 
RUNNING ANGLE ——— 
Dr) / 


BALANCE BETWEEN 
MODEL AND SPRING 


TYPICAL YAWING 
MOMENT CURVE 


FIG. 4 





SPRAY ROOT LINES 


Py 
= —— 


—e | 


FIG.LS 





REFERENCE Ra TANCt FOR 
AFTERBODY BOTTOM 


| | — yu 
| | 








2°94 


90+) 

ey2>"9 

NOLLIGNOD ONINYId Ni 

~ TINH NO INDY S3IN0s 


IWA0Nd IAUM 
“WNIQALISNOT 


921-2, 


—, 





“8 'br 


0568 


4 7 


00% th 
% 








ae 














6°94 
Guady dalla 


NV Id AGOGGd LAV cor oe 





692249 AWea%D hE 2,49 
a 3d 


<\ 4 


hE CG 2 AD | 
Qhes “D | 
69¢6=%) | 
Any 
| “ dais 
| 
| 
< e tT i 





AU40dd SAVM O3LNdWOD 


LHDIaH aAWM G3LSnrau a 





SE9= 4h O1°9l4 
egg: ND 
go: Y%9 








HULL FORCE DIAGRAMS 


FOREBODY AFTERBODY 
20,650 ecg 4475 ANTS 


8570 












Cy = 2.63 
21,250 
470 (a) 210 +220 
14,050 19,050 6175 6175 
Re Cy = aot 3436 
5,134 
19,678 
6459 
5,890 
(b) 
551 5402+ 290 
17,750 17,750 7425 [ 
Cy= 234+ 





514 FIG. 11 361+ 328 





YAWING MOMENT - INCH PCUNLS 


5.0 


4.0 


5.0 


2.0 


j.o 


SPRING CALIBRATION 


YAWING MCMENT VS AYAW 


RED SPRING 
5.54 INiBS/ DEGREE 


YELLOW SPRING 
2.414 INLBS/ DEGREE 


-2 4 .6 8 4.0 ie 


DB YAW #5 .E - DEGREES 
FiG.i2 





PORT 





YAWING MOMENT VS YAW ANGLE 
FREE-TO-TRIM 


Cv = 2.63 
Ca = 0.8 
Cm *¥O 


© 
© 
© 
o yh. 
oN 
——~- =-0038 
90 


-.O)} 


~.02 


-.03} 


~.04 


.QO3 


DISPLACED AXIS DUE TO 
ERROR IN YAW SETTING 


sé 
hp 44--™ 
© 
© 
© 
2 ty b 
a, =e 
~ 
NN Beco 
NY 
ae 
~ 


FIG.13 


STBD 





YAWING MOMENT VS YAWANGLE WV 
FIXED TRIM (12.6°) 


FRONT Pivots (NZ) 05 
Cv = 2.63 
Ca = 0.6 
OF 
03 © 
\ 5 
e- sy oe 
\ \ 
| f. \e O48 
| © © y \e 
| port STBD 
8 & yy. 2 NJ 2 4 6 8 
| , 
ON2 ~ ~.0075 —.04L> at © 
aA& 


mR 


-.0 oO \ 
aN 


AXIS DISPLACED DUE TO 
93). ERROR IN YAW SETTING 


FIG.14 








Q 
YAWING MOMENT VS YAW ANOLE. NV 
FIXED TRIM (12.6°) 


BACK PIVOTS (N,’) 6 
Cv = 2.63 
Ca « 0.8 

| O%F © 

03 


es - 


oO a © 
PORT STBD 


8 0 uF 5 N 2° if 6 e 
ON oOo @Mo 
t 
ON, zs - QOQ50 0% Q © © 
O® \ 
| AXIS SHIFT DUE TO 
ERROR IN YAW SETTING 

-.02 
-03 
-.0% 


FIG.15 





YAWING MOMENT VS YAW ANGLE 


FREE -TO-TRIM 
tv Cv = 2.46 
Cor 0.8 
Cu = O 


<) 


PORT 


‘ 
ON 2.0037 
oO -.08 


-02 


-.03 


03 


AXIS DISPLACED DUE TO 
ERROR IN YAW SETTING 


FIG.16 


16 


t% 


12 
© 410 
9 
STBD 
& 8 





YAWING MOMENT VS YAW ANGLE 


FIXED TRIM (Ceil.2°) 
BACK PIVOTS (Ni ) 


Cv = 2.46 
CorO8 
) 
© 
© 
eae 4 
. 
PORT 
b uf 2 


O¥ 
OS 
Oe 


.O4 


~.O1 


~.02 


» 


-03 


AXIS DISPLACED DUE TO 
ERROR IN YAW SETTING 


STBPD 


FIG. 17 





YAWING MOMENT VS YAW ANGLE 


FIXED TRIM (% = 11.27) 
FRONT PIVOTS (N32) 





Cv = 2.46 .O5 
Caz 0.8 
04 
P) .03 
- ry 
02 
9 
\ 
: ° .O1 
DUE Oe cRROR 7K eo \o 
PORT HEEL SETTING ___¢ STBO 
8 o Ue 2 \@ 2 4 b 8 
DNe WS 
2 =-,0052 aS 
d® -.01 ‘ 
AXIS DISPLACE MENT a \ 
DUE TO ERROR IN NS 
YAW SETTING : 
atOz sd 
-.93 
~.O4 
~.05 


FIG.18 





YAWING MOMENT VS YAW ANGLE 
FREE-TO- TRIM 


3 Cer 2.37 1o 
Co= 0.8 Cj 
Cu + O C) 






Me 
% 
N Cc 12 
~ 10 
Tt 
06 
\ 
o+ 
\e 
6 J 
SS sale 
\ () q 
PORT | STBD 
8 © uP ] 2 ub. b 2B 
\ 
Q © 
© \ 5 


—~02 


= - .O130 


St. 
gr 


«> 


-.08 \ 


FIG.19 





YAWING MOMENT VS YAW ANGLE 
FIXED TRIM (T= 10.5") 
BACK PIVOTS (N/) 








Cy = 2.34 
Co 70.8 
; 
N 
~ 
PORT 
8 © 4 
§ 
ON =—.Q130 
ae 
Oo 
©\o - 
© 


-.06 


FIG.20 





YAWING MOMENT VS YAW ANGLE 


FIKED TRIM (T <10.5°) 
FRONT PIVOTS (N2’) 


Cy = 2.54 
Ca *+O0.8 
N 
06 
SN 0% 
\ ‘a 
SS 
~. : 
; 02 
~ ; 
2 \ 
ny \ 
~ 
fs x 
5 2 Mie 7 a 2 3 
N 
1 
\ 
ON, = -.018 ae > | 
de® . . 
aN 
oO » 
- oF © 
ERROR IN 
YAW SETTING 
~ 0b 


FIG. 2i 





STATIC STABILITY DERIVATIVES 
VS SPEED COEFFICIENT 


FOR COMPARISON OF TESTED 
AND COMPUTED VALUES 


a Ov4e (IN RADIANS) 


i 
Ore (IN RADIANS) 


-0.5 
-O4] .2 
-0.3 


a 0.2 








d 
1 ON computen 
ae 


2.2 2.3 2.4 2.5 2.6 2.7 
SPEED COEFFICIENT ( Cv) 


FIG.22 





| Thesis 


17305 
| DIQ De Callies 


Evaluation Of static 
Stability derivatives 


From standard yawing 
Eosts 





Evaluation of static stability derivativ — 






